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I. INTRODUCTION 

Ultra-relativistic heavy-ion collision experiments pro¬ 
vide enough energy to cross the threshold where hadronic 
systems are converted into a novel state of matter, the 
quark-gluon plasma (QGP) [HE]- In this medium, par- 
tons acquire a much longer mean-free path and can move 
more freely due to asymptotic freedom mz]. Below the 
transition temperature, and up to a few times Tc, Quan¬ 
tum Chromodynamics (QCD) is essentially nonpertur- 
bative. For vanishing baryon chemical potential. Lattice 
QCD provides an excellent description of the full ther¬ 
modynamics (see e.g. Refs. 011] for recent results) that 
can be compared to improved perturbative calculations 
for temperatures a few times and higher. However, 
for nonzero baryon densities solid lattice results are still 
precluded by the well-known Sign Problem uni, and one 
has to resort to a combination of perturbative calcula¬ 
tions and effective models. 

Using perturbative techniques to calculate the ther¬ 
modynamic potential has a long tradition in QCD |11L 
114] . Since then, analytic computations of the thermody¬ 
namic potential were performed over the years (see e.g. 
[11 [T51 - ET] !. and are currently known up to 0{g^\ng) in 
the coupling g at high temperatures and small chemi¬ 
cal potentials pT1126| . It turned out that the pure weak 
coupling expansion converges badly. The next-to-leading 
order (g^) shows a surprisingly good agreement with lat¬ 
tice data whereas the next contribution (g^) even causes 
a flip in the sign of the pressure. It was suggested that 
the bad convergence is initiated by the contributions of 
soft momenta, Psoft ~ gT. Therefore the pressure of hot 
QCD is split into two parts pqcd = Phard+Psoft [27j, where 
Phard ~ ttT are the hard modes. The poor convergence 
can be assigned to Psoft, whereas phard exhibits a good 
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convergence. There are ways to reorganize the perturba¬ 
tive series with a special focus on the soft sector. One 
successful approach is Hard Thermal Loop perturbation 
theory (HTLpt), see Refs. [3SHS2]- Another method is 
dimensional reduction |33l [34], which proceeds by first 
integrating out the hard modes and leaves an effective 
three-dimensional theory, dubbed EQCD. A very recent 
investigation of the Equation of State of quark matter 
is based on both of the mentioned approaches and guar¬ 
antees access to all values of temperature and density 
at order g^ [3S]. Besides that, at order p®, intrinsic in¬ 
frared problems appear and can only be solved by non- 
perturbative methods |36j. 

Although the first calculations of the thermodynamic 
potential with massive quarks also date back to the 70s 
[MEZIISZIISE], mass effects on the pressure were regarded 
as being negligible for two decades. The issue was recon¬ 
sidered at T = 0 in Ref. [39], with a modern MS scheme 
description, and provided corrections that can reach 20%, 
being relevant for the physics of neutron stars. Later, the 
two-loop result of Ref. (35] was extended to order g'^ in 
Ref. [40], which represents the current state-of-the-art 
perturbative calculation at T = 0 with massive quarks. 
The case at zero baryon density and finite temperature 
was investigated in Ref. m with consequences to quark 
mass thresholds (see also Ref. [42]). Mass effects on the 
Yukawa theory were also investigated for massive quarks 
and mesons in Refs. [331145] . In all these cases, results 
exhibit deviations 20-30%. Quasi-particle models that 
include finite quark masses were also considered [46] . 

In this note, we present results for several thermody¬ 
namic quantities within the next-to-leading order calcu¬ 
lation (NLO) of the thermodynamic potential in pertur¬ 
bative QCD at finite temperature and chemical potential, 
including non-vanishing quark masses. We are interested 
in the investigation of effects of consistently introduced 
massive quarks, which were surprisingly less studied in 
the past. A model that includes all of these depen¬ 
dences in a reasonable way is the p^-corrected thermo¬ 
dynamic potential with mj J 0. As already mentioned, 
it is known for a long time that the pure weak coupling 
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expansion converges badly. However, it is also known 
that the NLO corrected pressure is relatively close to the 
lattice data which is of course only accidentally. Never¬ 
theless in terms of a good agreement with lattice simu¬ 
lations the NLO computation provides a baseline from 
which to study further corrections. What is not known 
nowadays is the impact of the finite bare quark masses 
on the thermodynamics. The influence of the strange 
quark mass for example to the thermodynamical quan¬ 
tities was ignored so far in the literature even though 
the strange-quark plays an important role, regarding the 
temperatures achieved in present ultra-relativistic heavy- 
ion collision experiments. The question is how good is the 
convergence of the pressure in terms of the quark mass 
dependence. In Ref. [H] it was shown that the conver¬ 
gence is better and the results should be trustable already 
at NLO. Furthermore, since quark mass effects in psoft 
only appear through the thermal masses, Debye-masses 
Tod, bare quark mass effects should be well described 
in NLO. The only other term that may have a signifi¬ 
cant quark mass dependence is the NNLO contribution 
in Phard, but sincc the convergence of Phard is better this 
contribution is small. Indeed, we will show in the follow¬ 
ing that the model behaves reasonably well in terms of 
mass and radiative corrections. In order to accomplish a 
consistent calculation that also includes a running cou¬ 
pling and strange-quark mass, the calculation at is 
done in the MS renormalization scheme. In principle the 
calculation of the thermodynamic potential was done by 
J. I. Kapusta in 1979, see Ref. m- However, it is based 
on an obsolete renormalization scheme and an explicit 
numerical evaluation of the thermodynamic potential is 
missing. 

We profit from previous calculations of the thermody¬ 
namic potential to study other properties such as the en¬ 
tropy, the speed of sound and trace anomaly with massive 
quarks. These results are compared to lattice data and 
to higher-order optimized perturbative calculations with 
the purpose to investigate the trend brought about by 
mass corrections. Partially included higher order effects 
that enter our framework by implementing the running 
of the coupling constant and the strange-quark mass re¬ 
sult in an improved description of lattice data already at 
NLO. 


The paper is organized as follows. In section H we sum¬ 
marize the well-established calculation of the thermody¬ 
namic potential to this order. In section HI we present 
our results and compare them to lattice data and HTLpt. 
Section IV contains our summary. Appendix A presents 
some technical details of the renormalization procedure 
for nonzero quark masses. 


II. THERMODYNAMIC POTENTIAL 


The next-to-leading order thermodynamic potential 
has the diagrammatic form [47] : 



+ diagrams with counterterms + 0(3 loops), 

( 1 ) 

where solid lines represent fermions, curly lines gluons 
and dashed lines ghosts. Notice that including nonzero 
quark masses, the fermionic contributions can not be 
rescaled by the flavor number Nf. 

The only diagram that will be affected in a nontriv¬ 
ial manner by nonzero quark masses is the exchange di¬ 
agram. To make the discussion more self contained, we 
present in some detail the derivation of this contribution. 
In the Feynman gauge m, we have 



= \pVNGg^ f 
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(27r)3 J 


dh r 

(27r)3 J 


(Pk 

(27r)3 


x{2Tr)^ 6 {p-q-k)T^ ^ 
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( 2 ) 


where Nq = - I is the number of gluons, p^ = {p° = 

+ ii,p) is the fermion four-momentum {q^ is defined 

analogously) and k^ = = iuj^^,k) is the gluon four- 

momentum, with = 2n7rT and = (2n + l)7rT, n 
being an integer. The mass of the respective quark fla¬ 
vor is 771/. The trace is performed in Dirac space. It 
is straightforward to show that Eq. © (see e.g. Refs. 
HH HH] for details) reads: 


d^pd^q 


J (27r)3 EpEqLUpq 
X ^J^UJpqYji + J_U)pq'^2 
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with the definitions 
J^ = Ng 


+ p - EpEq 


El-u:lq 

Nf{p) = nf{Ep + p) + nf{Ep - p), 

Nf{q) = nf{Eq + p) + nf{Eq - p), 

51 = nf{Ep + p)nf{Eq + p) +nf{Ep- p)nj{Eq - p), 

52 = nf{Ep + p)nf{Eq - p) +nf{Ep- p)nj{Eq + p), 

( 4 ) 


where E^ = Ep± Eq, ujpq = yj\p-qf, miu}) = (exp(/3a;) - 
1)“^ and nf{E ± p) = (exp(/?(ii^ ± p)) + 1)“^. 

As customary, vacuum contributions can be absorbed 
by a constant shift of the thermodynamic potential. The 
ultraviolet-divergent terms 
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(27r)6 EpEqUJpq 

[J+ {E.+UJpq) - J- {E+ -UJpq)]Nf{p), 


( 5 ) 


Lb = 


-PV9‘f 


d^p<Pq 
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(27r)6 EpEqUJp 


-2J.E+nb{ujpq), (6) 


are renormalized in the MS scheme to (see Appendix [a| 
and Refs. |44l |48] for details) 


-REN 
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{2tt) 3 Ep 
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= 0 , 


( 7 ) 

( 8 ) 


so that the renormalized exchange contribution has the 
form 

a r°° r ^ 

f^exch= 7 ^ / dEp / dEq / d{cos9) 

47 r'^ Jrrif Jmf J -1 

X yjE^-m'jyjE^ - TO^{J+Si + J_E2} 


NgasT^ 

Gtt 

Ncasm'j 


where 9 is the angle between the momenta p and q. Equa¬ 
tion recovers previous results [M M in their lim¬ 
its provided one corrects for the renormalization scheme. 
Note the explicit scheme dependence at this order due to 
the nonzero quark mass. The integrals can not be cast 
in analytic form as in Ref. [31] so that they have to be 
solved numerically. 

For the complete thermodynamic potential to this or¬ 
der, one has to add to the exchange contribution the 
following well-known terms 






dE{E^-m)fENf{p), (10) 

J m 


= -'—NgT^^ + —NNgT'^. 

glue G 36 = ^ 


III. RESULTS 


( 11 ) 


To investigate the influence of nonzero quark masses 
for the thermodynamics, we fix the up and down down- 
quark masses to 


= 2.3MeV and = 4.8MeV, 


( 12 ) 


and incorporated the running of the coupling and strange 
quark mass following Ref. [39], so that 


47r 
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(14) 


dEpJEj-m)Nf{p) 

Jmf V 

L, ‘iS,V'EJ-™5«,{p)[2 + 31n(Aj 


where L = 21n(A/A^), /3o = 11 - 2Nfl3 and /?i = 51 - 
19A^//3. The scale A^ and the invariant mass rhs are 
fixed by requiring 0.3 and rrig -22 100 MeV at A = 2 
GeV [^; one obtains Aj^ 22 380 MeV and rhs - 262 
MeV. With these conventions, the only freedom left is 
the choice of A. The renormalization scale is usually 
chosen to be A = 2\J 4- p^. The band uncertainties 
in our plots shown below result from the variation of A 
by a factor of two. In all calculations the quark chemical 
potentials are set to be equal pf = p = psl^- In order to 
cross check our computations the results are compared to 
the case of zero temperature and zero chemical potential. 


A. Pressure, energy density, entropy and trace 
anomaly 

We start with the zero-temperature case. The one mas¬ 
sive flavor case can be computed analytically, yielding, in 
the MS scheme: 


( 9 ) 
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(15) 


where u = 
Figure 


1 ^ 


p^ - mq. The pressure is given by p = -fl. 
shows a comparison of calculations for the 
pressure at different orders for 2-1-1 flavors at vanish¬ 
ing temperature. As customary, we plot the Stefan- 
Boltzmann-normalized (SB) pressure. The red band is 
taken from Ref. [ID] (see also Ref. [SO]), where the au¬ 
thors assumed massless light quarks but a non-vanishing 
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2+1 and 3 Flavors, T=0 MeV 



Mb [GeV] 


2 Flavors, T=300 MeV 



FIG. 1. Comparison of the normalized pressure for differ¬ 
ent orders in perturbation theory; m/ + 0 means that quark 
masses were chosen according to Eq. (121 and Eq. (14l. The 
red band is taken from Ref. m (see also Ref. [SS]) 


strange quark mass (The light quark mass effects are, of 
course, expected to by minor.). This plot quantifies the 
effect from higher order contributions. As expected, at 
vanishing temperature they are less severe, even though 
the pressure is non-ideal even for very high densities. 


2+1 and 3 Flavors, T=300 MeV 



In Figure we show the normalized pressure for dif¬ 
ferent numbers of flavors varying the baryon chemical 
potential and keeping the temperature fixed at T = 300 
MeV. We illustrate the results for different quark masses: 
the upper plot in Fig. shows the difference in the nor¬ 
malized pressure pjpsB between the results for massless 
quarks and those for non-vanishing quark masses with 
mf = 200 MeV (a commonly encountered mass scale in 
lattice calculations). The difference in the normalized 
pressure is roughly - 9%, so that direct comparisons to 
lattice data with unphysical quark masses should be done 
with caution. The lower plot in Fig. shows the differ¬ 
ence in the normalized pressure pjpsB between the mass¬ 
less case and the physical mass case for 2-1-1 and 3 flavors. 
This difference in the normalized pressure is less than 1%, 
which shows that temperature effects easily dilute mass 
effects that are stronger in the cold and dense case (see 
Refs. [39l HD])- Both of the plots in Fig. unfold the 
behaviour of the pressure for a number of 2 (top) and 3 
(bottom) flavors with all have the same mass of 1 GeV. 


FIG. 2. Normalized pressure for 2 (top) and 2+1 / 3 flavors 
(bottom) for different quark masses; m/ ^ 0 means that the 
quark masses were chosen according to Eq. ( |12[ ) and Eq. ( |14[ |. 

MeV and ps = 400 MeV. In both plots a drastic devia¬ 
tion again appears for increased quark masses of to/ = 1 
GeV. 


In Fig.|^the normalized pressure for pB = 0 (top) and 
for pb = 400 MeV (bottom) is depicted and compared 
to lattice simulations m and results of other perturba¬ 
tive investigations that include higher-order effects, for 
illustration. Of course, higher-order corrections are ex¬ 
pected to modify appreciably our simplified description. 
However, taken as a toy model, this description is not far 
from lattice results within the interesting regime in tem¬ 
perature. The green bands on the plots of Fig. again 
illustrates the large effects that can be brought by the 
inclusion of heavier quarks. 


We investigate the same situation where the normal¬ 
ized pressure is plotted over the temperature in Fig. 
This time the difference for the 2-flavor case with in¬ 
creased masses of to/ = 200 MeV (upper plot in Fig. [3) is 
larger. For the lowest illustrated temperature of T = 150 
MeV it amounts to almost ~ 70% at very low tempera¬ 
ture. Note that at this temperature and for pB = 400 
MeV as ~ 0.54. For 3 flavors with physical quark masses 
compared to massless quarks (lower plot in Fig. there 
is also an increase that is approximately 11% at T = 150 


'inally, another comparison can be made by introduc- 
the pressure difference 


Ap = p{T, A, p) - p{T, A, 0). 


(16) 


In Fig. we show the so-called scaled pressure difference 
which is the pressure difference, Eq. (16), normalized by 
T^, which eliminates all contributions to the pressure 
that do not depend on p. Notice that lattice data only 
includes corrections up to order p^. Again, taken as a 
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2 Flavors, 1^8=400 MeV 


2+1 and 3 Flavors, Hb= 300/400 MeV 




3 Flavors, iib= 400 MeV 



FIG. 3. Normalized pressure for 2 (top) and 2+1 / 3 flavors 
(bottom) for different quark masses; m/ + 0 means that the 
quark masses were chosen according to Eq. (|12[) and Eq. (|14[). 


2+1 and 3 Flavors, |Jb= 0 MeV 



2+1 and 3 Flavors, iib= 400 MeV 



FIG. 4. Comparison of the normalized pressure for 2+1 / 3 
flavors, hb = 0 (top) and fiB = 400 MeV (bottom) with HTLpt 
|32 |. dimensional reduction [5S], and lattice results, denoted 
by WB, [511152j : nif t 0 means that the quark masses were 
chosen according to Eq. (|12[) and Eq. (|14[|. 


FIG. 5. Comparison of the scaled pressure difference 
with HTLpt [32] and lattice results, denoted by WB, | 51| : 
ruf t 0 means that the quark masses were chosen accord¬ 
ing to Eq. ( |T2| ) and Eq. |T4| ). 


simplified model, our description captures well the de¬ 
pendence on /r including quark mass effects of this quan¬ 
tity as a function of temperature. 


One can calculate the energy density from the expres¬ 
sion of the pressure potential using the relation 


dp dp 
dT dp ^ 


(17) 


Note that we fix A first and then take the derivative, 
in contrast to e.g. Ref. |31], where the derivatives are 
taken first and A is fixed afterwards. In Fig. we plot 
the normalized energy density versus the temperature for 
vanishing baryon chemical potential (top) and for pB = 
400 MeV (bottom). The interpretation of the results is 
analogous to the one for the pressure. Here we also show 
the effect from a heavier quark to the energy density. 
Finally, the entropy can be directly computed from the 
pressure, S = dpjdT, and is exhibited in Fig.j^for ps = 0 
(top) and pb = 400 MeV (bottom). 


The trace anomaly of QCD [29] is defined as X = e-3p. 
It is related to deconfinement and the gluon condensate 
[53] ■ therefore non-perturbative in its origin. This quan¬ 
tity vanishes for a massless, noninteracting gas since, in 
this case, e = 3p. The mass dependence of the results 
for the trace anomaly could be of special interest since 
nonzero masses break explicitly the scale invariance of 
QCD. On the other hand, interactions also break scale 
invariance as can be seen in the running of as- We com¬ 
pute X for 2+1 and 3 flavors, and illustrate it in Fig. for 
Pb = 0 (top) and pb = 400 MeV (bottom). Our results 
systematically underestimate the trace anomaly, show¬ 
ing that higher order terms are necessary to capture the 
full running effect. We show also results for a higher 
quark mass to quantify how it affects this quantity. 
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2+1 and 3 Flavors, |Jb= 0 MeV 



2+1 and 3 Flavors, 


2+1 and 3 Flavors, |Jb= 0 MeV 





FIG. 6. Comparison of the normalized energy density for 
fj-B = 0 (top) and /is = 400 MeV (bottom) with HTLpt [32] 
and lattice results, denoted by WB, [nHSH; rrif t 0 means 
that the quark masses were chosen according to Eq. (12 I and 

Eq. pi. 


The relative contributions of the mass and the running 
Qfs to I are illustrated in Fig. It is obvious that the 
one from running clearly dominates. The result for 
a fixed but a nonvanishing strange quark mass of 95 
MeV is rather small. The sum of both contributions is 
also plotted. 


B. Speed of sound and susceptibilities 


The speed of sound Cg is given by 

2 dp 


(18) 


It is a measure of the stiffness of the equation of state, 
and is naturally bounded due to causality Ug < 1 and ther¬ 
modynamic stability > 0. It is relevant at low densities 
and high temperatures in the hydrodynamical descrip¬ 
tion of the evolution of the quark-gluon plasma formed 
in relativistic heavy-ion collisions |54l - f56| . On the other 
hand, the way the squared speed of sound approaches 
asymptotically 1/3 is relevant in the astrophysical con¬ 
text, as discussed in Ref. since there are scenarios 


FIG. 7. Comparison of the normalized entropy for fj,B = 0 
(top) and hb = 400 MeV (bottom) with HTLpt [32] and lat¬ 
tice results, denoted by WB, [SmSD; mf t 0 means that the 
quark masses were chosen according to Eq. (121 and Eq. (141. 


2+1 and 3 Flavors, Pb= 0 MeV 




T[GeV] 


FIG. 8. Comparison of the trace anomaly for hb = 0 (top) and 
Hb = 400 MeV (bottom) with HTLpt |32| and lattice results, 
denoted by WB, | 511I52 |: m/ + 0 means that the quark masses 
were chosen according to Eq. (12 I and Eq. (14l. 
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2+1 / 3 Flavors, Hb= 400 MeV 2+1 and 3 Flavors, Hb= 0 MeV 




FIG. 9. Comparison of the trace anomaly contributions for 
Hb = 400 MeV 


2+1 and 3 Flavors, T=0 MeV 



Fb [GeV] 


FIG. 10. Speed of sound for 2+1 / 3 (massive) flavors at 
T = Q\ ruf t Q means that the quar k masses were chosen 
according to Eq. (121 and Eq. (14l. 


within neutron stars which seem to violate < 1/3 (see 
also Ref. [55|). 

The speed of sound at zero temperature is shown in 
Fig.[^ which shows that physical quark mass effects and 
a running coupling constant do not result in a sizable shift 
from the limiting value = 1/3. In Fig. 11 we show the 


computation of for non-vanishing temperature at /is = 
0 (top) and fiB = 400 MeV (bottom). Again, the effect 
from nonzero masses become relevant only for heavier 
quarks. The band that is caused by the variation of A by 
a factor of two, is only plotted down to densities of /is ~ 
2, 7 GeV. Below that density the lower line of the band 
starts to increase rapidly, which is of course no physical 
effect, because as does so due to the logarithm in L = 
21n(A/A^), see Eq. (131. 


2+1 and 3 Flavors, |Jb=400 MeV 



FIG. 11. Comparison of the speed of sound for fiB = 0 (top) 
and fiB = 400 MeV (bottom) with HTLpt |32], and lattice 
results, denoted by WB, iH |S2]; nif 0 means that the 
quark masses were chosen according to Eq. (121 and Eq. (141. 


In Fig. 
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_ iCP 
Vuuuu ~ ' 


For instance, xt - Xv 

we plot xt normalized to its SB-limit and 


compare it with several other data. Additionally we plot 
the case for 3 flavors with equal masses of 1 GeV. In 
Fig. 13 the second- (top) and the fourth-order (bottom) 
baryon number susceptibilities, given by 


XUT) 


d^p 


dp- 




( 20 ) 


are displayed. Whereas for physical quark masses the 
mass effect on the susceptibilities is arguably small, for 
heavier quarks the effect is dramatic, even at comparable 
and large values of the temperature. Since these quan¬ 
tities seem to be much more sensitive to mass effects, 
comparisons to lattice results with unphysical high quark 
masses should be made with caution. 


We also compute quark number susceptibilities, which 
are defined as derivatives with respect to the correspond¬ 
ing quark chemical potentials jl = {pu,pd, ■ ■ ■, pNf ) as 


... 


( 19 ) 


IV. SUMMARY 

We presented a systematic study of the effects from 
nonzero quark masses on the thermodynamics of pertur¬ 
bative QCD to next-to-leading order at non-vanishing /i 








































1 Flavor, Hb= 0 MeV 



(unphysical) quark masses has a limited validity mean¬ 
ing that the impact of quark masses is not negligible. 
In this vein, susceptibilities seem to be particularly sen¬ 
sitive to higher quark masses. Although a calculation 
at higher order (including resummations) is imperative 
for any quantitative statement, our results point to the 
role heavier quarks will play when incorporated as in the 
case of the early universe where charm quarks appear at 
a few times T^. The interplay between quark masses and 
isospin is also a relevant issue (see Ref. [60]) and should 
be addressed in forthcoming work. 


FIG. 12. Comparison of the normalized 4th order u- 
quark susceptibility with HTLpt [32] , dimensional reduc¬ 
tion m and lattice results, denoted by WB, m and 
BNL, |59| : m[ t Q means that the quark masses were 
chosen according to Eq. (121 and Eq. (141. 


2+1 and 3 Flavors, |Jb= 0 MeV 



2+1 and 3 Flavors, |Jb= 0 MeV 



T [GeV] 


FIG. 13. Comparison of the second (top) and fourth (bot¬ 
tom) order baryon number susceptibility with HTLpt, | 32 |. 
dimensional reduction, IS], and lattice results, denoted by 
WB, [52]; ruf t Q means that the bare quark masses were 
chosen according to Eq. (121 and Eq. (14l. 
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Appendix A: Renormalization 


We discuss the renormalization of L/ and Lb, for de¬ 
tails see Ref. [H] . The expressions derived in ([^ and ([^ 
are 


Lf=-pVg^ J 


dPpdPq 1 

(27r)6 EpEqUJpq 


X ■*" ^pq) J- (-^+ ^pg)] 


(Al) 


Lb = -pvg^ f - 2j.E,n,{u;pq). (A2) 

As described in Ref. [33] we define the auxiliary functions 

\A f 4,, _ dq'^dk"^ 2m^-p^qp 

(2^)2 ((g4)2 + £;2)((^4)2+^2^) (A3) 

X 2TrS{p‘^ - - k^), 


and T. We investigated the pressure, scaled pressure dif¬ 
ference, energy density, entropy, trace anomaly, speed of 
sound and susceptibilities and compared our findings to 
Hard Thermal Loop and Dimensional Reduction calcula¬ 
tions at vanishing quark masses and recent lattice data 
at physical quark masses. 

A comparison between perturbative results with van¬ 
ishing quark masses and lattice results with increased 


M (k'^')= 2m'^-p^qp, 

{27r)^ {{p^y + E^){{q^r + E^q) (A4) 

X 27ri5(p'^ - g^ - ^^), 

with p'" = (p°,p) = {ip‘^,p) and g'" = {q°,q) = 

The flavor index of the quark mass m/ is supressed in 
this appendix. Plugging in the integral representation 
of the Dirac delta function we factorize and solve the 
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integrations independently m- 

The function M. f{p^) can be written as 


Mf{p^) = 


AEqOJpq 


and as 


2m^ + ip^Eq + pq 2m^ - ip^Eq + pq 

ip^ + Eq + iOpq —ip^ + Eq + iOpq 

(AS) 




AEpEq 


2E+ 


El - {ik^y 


(A6) 


Eq. and Eq. ( |A2[ ) can be expressed in terms of the 
auxiliary functions as 


Lf = pVg^ f 

4.Ng 




E„ 


(27r)9 

Nf{p)Mf{-iEp), 


(A7) 


Lb = 






{2nf 


ANcribiuj) . , . . ^ 


(A8) 


Using the expression for M.f{p'^) and Mb{k'^) one arrives 
at 

-Wo^iVr f f f 

^ ^ ^ J {2ttY \J f27r')6 J 


Lf = 


X (27r)^(5(p- q- k)2'iT5{p'^ - - k^) 




((g4)2+K2)((fc4)2+^2) 


p‘^=-iEp 


(A9) 


U-msNaJ 7^—\J j^J 7^ 


(27r)3 uj (-/ (27r)9 J (Stt)^ 

X {2'KyS{p- q- k)2'K6{p'^ - - k'^) 

217^ -p^qp 

^(p4)2+R2)((q4)2+K2) 


k^=-iuj 


With the help of the definitions of and = (k^,k) = 
{ik‘^,k) it follows that 


X (-l){2TT)y6^^\p^ - q^ - k^^) 

2ny - pt^q^ 


(All) 


((g4)2+£;2)((fc4)2+^2) 


P=Ep 


Lb = -ApVg^NG f 


drk nb(u!) ( f drpdrq 


I 


(27r)3 w [J (27r)^ 

X (-l)(27r)4i(5l4)^pM - qt^ - fc^) 

277^-p^^qp ] 

^(p4)2+p;2)((^4)2+^2)|^^^^- 


(A12) 


The expressions in the curly brackets are just the ampu¬ 
tated (AMP) self-energy diagrams of a quark or a gluon, 
respectively, in the vacuum (VAC) where their momenta 
are put on-shell (M.S.): 




VAC 


'AMP 

M.S. 


= \-9 


/ 


d'^qtPk 


(27r)8 

X (27r)V‘*)(p^-g'^-fc^) 
^ Ng 2m^-p^g^ 

Nc m{q^ - m?-)kP- 


(A13) 


( \ VAC 

j = |-5 


/ 




(27r)8 

X (27r)4d(^)(p^-g^-fc'^) 

AT AT ~ 

X NoNf 


(p2 _ m‘2^(^q2 _ ^2^ 

(A14) 


fc2=0 


Defining p^^^ e (iw^^ + p,p) and q^^^ = (iuj^^^q) with 

Nf{p)-l 


r^Tr 


- m 


= 4to- 


2Eri 


ty 


rin ^new 


2/7^ + 1 
2a; ’ 


(A15) 


the vacuum self-energies can be written as 


f (-l)Tr 

Knew 


=9‘ 


/ 


d^p 

(27r)3 


Nf{p) - 1 
2E„ 




. VAC n 


'AMP 

M.S. -■ 


(AlO) xjy' ^^(27r)V^)(p^-g'^-fc'^) 


Ng 8m^ -4p^g^ 

t2^^2 


Ne (g" 


p=r] 


(A16) 


lUnew dnew y >/ /AMP 

M.S. 


=5 


r d^k 

prifa + n 

J (27r)3 

1. 2w J 




(A17) 
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8to^ - 4p^g^ 
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Additionally, one can transform it into the form: 


Lf = -2/3VNJi^ (- 1 ) 
I "^Pnew 


X Tr 




. VAC n 


f AMP 

M.S. 


MAT 


(A18) 


r 


4 

[ -^new <Znew \ V J J AMP 


AMP j A/fAT 
M.S. 


(A19) 


Lb 



(A21) 


The dashed box within the right diagrams indicates that 
this part of the diagram is put on-shell. Hence the UV- 
divergent part is completely covered by the dashed box. 
The matter part of the remaining integral does not con¬ 
tain any new divergences. 

Therefore the renormalization procedure reduces to the 
renormalization of the one-loop diagram in the vacuum. 
In the MS scheme the final renormalized expression for 
the amputated quark self-energy is: 




, VAC 


I AMP 
M.S.,REN 


(A22) 


img^Na 
(47r)2Ac Jo 


r 

Jo 


dx ^ In 


)( 2 - 


x) + {x- 1) 


whereupon the MAT index denotes that we take only 
the matter part of the expressions in the curly brackets. 
This means that the pure vacuum part has already been 
subtracted and, while and are evaluated on the 
mass shell, are not. These expressions can 

be assigned to the following two two-loop-diagrams 


where A = m^(l-x)^ and A is the renormalization scale. 
Finally one finds that 


rREN 


pVNcg^mJ 

dTT^ 


r dJp Nf{p) 

2 + 31n(- 

-]] 

J (27r)3 Ep 

L \r 

n)\ 


(A23) 

The renormalization of Lb poceeds analogeously but 
turns out to vanish. The general Lorentz structure of 
a bosonic self-energy in the vacuum reads 



= (k‘^gfj,^-k^k^)ll{k'^). 


(A24) 


In Eq. (A19) the bosonic self-energy is multiplied hygai^ 
which contracts with the Lorentz structure in Eq. (A241. 


(A20) The momentum of the ingoing gluon must be put on- 
shell, hence for a massless gluon k^ = 0. Subsequently, 
the whole expression vanishes because it is proportional 
to 


[1] E. V. Shuryak, “Theory of Hadronic Plasma,” 
Sov.Phys.JETP, vol. 47, pp. 212-219, 1978. 

[2] E. V. Shuryak, “Quark-Gluon Plasma and Hadronic Pro¬ 
duction of Leptons, Photons and Psions,” Phys.Lett., 
vol. B78, p. 150, 1978. 

[3] D. J. Gross and F. Wilczek, “Ultraviolet Behavior of 
Nonabelian Gauge Theories,” Phys.Rev.Lett., vol. 30, 
pp. 1343-1346, 1973. 

[4] D. Gross and F. Wilczek, “Asymptotically Free Gauge 
Theories. 1,” Phys.Rev., vol. D8, pp. 3633-3652, 1973. 

[5] D. Gross and F. Wilczek, “Asymptotically Free Gauge 
Theories. 2.,” Phys.Rev., vol. D9, pp. 980-993, 1974. 


[6] H. D. Politzer, “Reliable Perturbative Results for Strong 
Interactions?,” Phys.Rev.Lett., vol. 30, pp. 1346-1349, 

1973. 

[7] H. D. Politzer, “Asymptotic Freedom: An Approach to 
Strong Interactions,” Phys.Rept., vol. 14, pp. 129-180, 

1974. 

[8] S. Borsanyi, S. Durr, Z. Fodor, C. Holbling, S. D. Katz, 
S. Krieg, D. Nogradi, K. K. Szabo, B. G. Toth, and 
N. Trombitas, “QGD thermodynamics with continuum 
extrapolated Wilson fermions H,” Phys. Rev., vol. D92, 
no. 1, p. 014505, 2015. 




















11 


[9] H.-T. Ding, F. Karsch, and S. Mukherjee, “Thermody¬ 
namics of strong-interaction matter from Lattice QCD,” 
2015. 

[10] P. de Forcrand, “Simulating QCD at finite density,” PoS, 
vol. LAT2009, p. 010, 2009. 

[11] B. A. Freedman and L. D. McLerran, “Fermions and 
Gauge Vector Mesons at Finite Temperature and Den¬ 
sity. 1. Formal Techniques,” Phys.Rev., vol. D16, p. 1130, 

1977. 

[12] B. A. Freedman and L. D. McLerran, “Fermions and 
Gauge Vector Mesons at Finite Temperature and Den¬ 
sity. 2. The Ground State Energy of a Relativistic elec¬ 
tron Gas,” Phys.Rev., vol. D16, p. 1147, 1977. 

[13] B. A. Freedman and L. D. McLerran, “Fermions and 
Gauge Vector Mesons at Finite Temperature and Den¬ 
sity. 3. The Ground State Energy of a Relativistic Quark 
Gas,” Phys.Rev., vol. D16, p. 1169, 1977. 

[14] V. Baluni, “Nonabelian Gauge Theories of Fermi Sys¬ 
tems: Chromotheory of Highly Condensed Matter,” 
Phys.Rev., vol. D17, p. 2092, 1978. 

[15] S. Chin, “Transition to Hot Quark Matter in Relativistic 
Heavy Ion Collision,” Phys.Lett., vol. B78, pp. 552-555, 

1978. 

[16] J. 1. Kapusta, “Quantum Chromodynamics at High Tem¬ 
perature,” Nucl.Phys., vol. B148, pp. 461-498, 1979. 

[17] T. Toimela, “Perturbative QED and QCD at Finite 
Temperatures and Densities,” Int.J.Theor.Phys., vol. 24, 
p. 901, 1985. 

[18] P. B. Arnold and C.-X. Zhai, “The Three loop free energy 
for pure gauge QCD,” Phys.Rev., vol. D50, pp. 7603- 
7623, 1994. 

[19] P. B. Arnold and C.-x. Zhai, “The Three loop free energy 
for high temperature QED and QCD with fermions,” 
Phys.Rev., vol. D51, pp. 1906-1918, 1995. 

[20] C.-x. Zhai and B. M. Hastening, “The Free energy of hot 
gauge theories with fermions through g**5,” Phys.Rev., 
vol. D52, pp. 7232-7246, 1995. 

[21] K. Kajantie, M. Laine, K. Rummukainen, and 
Y. Schroder, “The Pressure of hot QCD up to g6 
ln(l/g),” Phys.Rev., vol. D67, p. 105008, 2003. 

[22] A. Vuorinen, “Quark number susceptibilities of hot QCD 
up to g**6 In g,” Phys.Rev., vol. D67, p. 074032, 2003. 

[23] A. Vuorinen, “The Pressure of QCD at finite temper¬ 
atures and chemical potentials,” Phys.Rev., vol. D68, 
p. 054017, 2003. 

[24] A. Ipp and A. Rebhan, “Thermodynamics of large N(f) 
QCD at finite chemical potential,” JHEP, vol. 0306, 
p. 032, 2003. 

[25] A. Ipp, A. Rebhan, and A. Vuorinen, “Perturbative QCD 
at nonzero chemical potential: Comparison with the 
large Nff) limit and apparent convergence,” Phys.Rev., 
vol. D69, p. 077901, 2004. 

[26] A. Ipp, K. Kajantie, A. Rebhan, and A. Vuorinen, 
“The Pressure of deconfined QCD for all temperatures 
and quark chemical potentials,” Phys.Rev., vol. D74, 
p. 045016, 2006. 

[27] E. Braaten and A. Nieto, “Free energy of QCD at high 
temperature,” Phys.Rev., vol. D53, pp. 3421-3437, 1996. 

[28] J. O. Andersen, M. Strickland, and N. Su, “Three- 
loop HTL Free Energy for QED,” Phys.Rev., vol. D80, 
p. 085015, 2009. 

[29] J. O. Andersen, L. E. Leganger, M. Strickland, and 
N. Su, “The QCD trace anomaly,” Phys.Rev., vol. D84, 
p. 087703, 2011. 


[30] J. O. Andersen, L. E. Leganger, M. Strickland, and 
N. Su, “Three-loop HTL QCD thermodynamics,” .JHEP, 
vol. 08, p. 053, 2011. 

[31] S. Mogliacci, J. O. Andersen, M. Strickland, N. Su, and 
A. Vuorinen, “Equation of State of hot and dense QCD: 
Resummed perturbation theory confronts lattice data,” 
JHEP, vol. 1312, p. 055, 2013. 

[32] N. Haque, A. Bandyopadhyay, J. O. Andersen, M. G. 
Mustafa, M. Strickland, et al, “Three-loop HTLpt ther¬ 
modynamics at finite temperature and chemical poten¬ 
tial,” JHEP, vol. 1405, p. 027, 2014. 

[33] T. Appelquist and R. D. Pisarski, “High-Temperature 
Yang-Mills Theories and Three-Dimensional Quantum 
Chromodynamics,” Phys. Rev., vol. D23, p. 2305, 1981. 

[34] P. H. Ginsparg, “First Order and Second Order Phase 
Transitions in Gauge Theories at Finite Temperature,” 
Nucl. Phys., vol. B170, pp. 388-408, 1980. 

[35] A. Kurkela and A. Vuorinen, “Gool quark matter,” 2016. 

[36] A. D. Linde, “Infrared Problem in Thermodynamics of 
the Yang-Mills Gas,” Phys.Lett., vol. B96, p. 289, 1980. 

[37] B. Freedman and L. D. McLerran, “Quark Star Phe¬ 
nomenology,” Phys. Rev., vol. D17, p. 1109, 1978. 

[38] E. Earhi and R. Jaffe, “Strange Matter,” Phys.Rev., 
vol. D30, p. 2379, 1984. 

[39] E. S. Fraga and P. Romatschke, “The Role of quark 
mass in cold and dense perturbative QCD,” Phys.Rev., 
vol. D71, p. 105014, 2005. 

[40] A. Kurkela, P. Romatschke, and A. Vuorinen, “Cold 
Quark Matter,” Phys.Rev., vol. D81, p. 105021, 2010. 

[41] M. Laine and Y. Schroder, “Quark mass thresholds in 
QCD thermodynamics,” Phys.Rev., vol. D73, p. 085009, 
2006. 

[42] X. Wang and J.-r. Li, “Effective two loop thermodynamic 
potential with fermions in the real time formalism of ther¬ 
mal field theory,” Commun. Theor.Phys., vol. 33, pp. 253- 
262, 2000. 

[43] L. Palhares and E. Fraga, “Cold and dense perturbative 
Yukawa theory with massive fermions,” Braz.J.Phys., 
vol. 37, pp. 26-29, 2007. 

[44] L. F. Palhares and E. S. Fraga, “Perturbative Yukawa 
theory at finite density: The Role of masses and renor¬ 
malization group flow at two loops,” Phys.Rev., vol. D78, 
p. 025013, 2008. 

[45] E. S. Fraga, L. F. Palhares, and M. B. Pinto, “Nonpertur- 
bative Yukawa theory at finite density and temperature,” 
Phys.Rev., vol. D79, p. 065026, 2009. 

[46] A. Peshier, B. Kampfer, and G. Soff, “The Equation 
of state of deconfined matter at finite chemical poten¬ 
tial in a quasiparticle description,” Phys. Rev., vol. G61, 
p. 045203, 2000. 

[47] J. Kapusta and C. Gale, “Finite-temperature field the¬ 
ory: Principles and applications,” 2006. 

[48] L. Palhares, “Materia densa e fria: Termodinamica, 
transicoes de fase e renormalizacao,” Dissertacao de 
Mestrado, (Universidade Federal do Rio de Janeiro, 
2008). 

[49] S. Eidelman et al, “Review of particle physics. Particle 
Data Group,” Phys.Lett., vol. B592, pp. 1-1109, 2004. 

[50] E. S. Fraga, A. Kurkela, and A. Vuorinen, “Interacting 
quark matter equation of state for compact stars,” As- 
trophys.J., vol. 781, no. 2, p. L25, 2014. 

[51] S. Borsanyi, G. Endrodi, Z. Fodor, S. Katz, S. Krieg, 
et al, “QCD equation of state at nonzero chemical po¬ 
tential: continuum results with physical quark masses at 



12 


order JHEP, vol. 1208, p. 053, 2012. 

[52] S. Borsanyi, G. Endrodi, Z. Fodor, A. Jakovac, S. D. 
Katz, et al., “The QCD equation of state with dynamical 
quarks,” JHEP, vol. 1011, p. 077, 2010. 

[53] E. Megias, E. Ruiz Arriola, and L. Salcedo, “Trace 
Anomaly, Thermal Power Corrections and Dimension 
Two condensates in the deconfined phase,” Phys.Rev., 
vol. D80, p. 056005, 2009. 

[54] S. Jeon and U. Heinz, “Introduction to Hydrodynamics,” 
2015. 

[55] C. Gale, S. Jeon, and B. Schenke, “Hydrodynamic Mod¬ 
eling of Heavy-Ion Collisions,” Int.J.Mod.Phys., vol. A28, 
p. 1340011, 2013. 

[56] U. Heinz and R. Snellings, “Collective flow and 
viscosity in relativistic heavy-ion collisions,” 
Ann.Rev.Nucl.Part.Sci., vol. 63, pp. 123-151, 2013. 


[57] P. F. Bedaque and A. W. Steiner, “Sound velocity bound 
and neutron stars,” Phys.Rev.Lett., vol. 114, p. 031103, 
2015. 

[58] A. Gherman, T. D. Cohen, and A. Nellore, “A Bound 
on the speed of sound from holography,” Phys.Rev., 
vol. D80, p. 066003, 2009. 

[59] A. Bazavov, H. T. Ding, P. Hegde, F. Karsch, C. Miao, 
S. Mukherjee, P. Petreczky, C. Schmidt, and A. Velytsky, 
“Quark number susceptibilities at high temperatures,” 
Phys. Rev., vol. D88, no. 9, p. 094021, 2013. 

[60] E. S. Eraga, L. F. Palhares, and C. Villavicencio, “Quark 
mass and isospin dependence of the deconfining critical 
temperature,” Phys. Rev., vol. D79, p. 014021, 2009. 

[61] P.Romatschke, “Thermal corrections to the sunset di¬ 
agram,” Diplomarbeit, (Technische Universitaet Wien, 
2001 ). 



